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We study the characteristics of the QCD phase transi- 
tions in dense hadronic matter using the Skyrme model con- 
structed on S 3 . We find numerically the localized solutions on 

c3 ■ S 3 using the extended Skyrme model which implements correctly 

the scale symmetry of QCD. The transition from the localized 
phase to the delocalized phase is found to be of first order at the 
critical radius of the hypersphere, L c . The chiral restoration and 
the gluon decondensation also take place at the same critical size. 
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The Skyrme model [Q describing baryons as topological solitons (skyrmion) 
arising from the effective meson theory of QCD in the large iV c -limit has been 
proved to be quite successful in calculating the various physical properties of 
baryons which are in fair quantitative agreements with experimental data. 
Given the success of the Skyrme model, there have been also interesting de- 
velopments to investigate the behavior of hadronic matter at higher density. 
One of the interesting methods, pioneered by Manton and Ruback^, ^| is 
to put one skyrmion on a hypersphere S 3 (p) with three dimensional volume 
2ii 2 p 3 . Although the obvious disadvantage of this method is that the physical 
space is not a hypersphere with curvature, this approach has been adopted to 
investigate the problems of dense matter in the Skyrme model motivated by 
the observation that it gives, in a simple frame work, similar results as in the 
lattice skyrmion approach Q and also by its mathematical simplicity. The 
hypersphere approach has been studied not only in the Skyrme model 0. [| 
but also recently in the Nambu-Jona-Lasino model 0. 

The skyrmion on the hypersphere of infinitely large radius is well local- 
ized so that it has the same properties as an isolated skyrmion in the flat 
space. The finite size solution which is well localized is possible due to the 
presence of the interaction term(Skyrme term or higher derivative terms), 
which prevents it from being collapsed. However, by shrinking the radius of 
the hypersphere, the effect of interaction term becomes much stronger than 
the kinetic term because of the curvature of the hypersphere. Therefore, it 
is possible that for the radius smaller than a certain critical radius p c , the 
localization mechanism due to the kinetic term becomes inactive and the 
skyrmion will spread out over the whole hypersphere uniformly. Manton || 
interpreted it as a signature for the quark-hadron phase transition. More- 
over, in relation to the chiral symmetry, the observation of vanishing order 
parameter (cr) at the critical density leads to the interpretation that it would 
be a signature for chiral symmetry restoration. 

In QCD, one of the important order parameter, is the gluon condensate 
(Ol-F^F^lO), which breaks the scale invariance in the hadron phase of QCD. 
The simplest way to implement this scale symmetry of QCD is to introduce 
a scalar field, xH> a dilaton or a glueball field in the effective Lagrangian. 
The scale anomaly of QCD in low energy physics is provided by the effective 
potential of x of the following generic form, 

fx 4 (lnx 4 -l), (1) 



where B is related to the bag constant. Eq.(|l]) has its minimum at nonzero 
Xo = (x)o which correctly reproduces the scale anomaly ||, and it is related 
to the gluon condensate, 



M 



Xo = ~ ^°I F ^|0> (2) 

Recently the scale symmetry property of QCD in the effective Lagrangian 
has been proved to be very useful in discussing the QCD phase transition JT0 



at finite temperature and finite density and also in investigating the medium 



effect on the parameters in the effective Lagrangian ||f f||. Gomm, et al[fl2 
studied the bag formation in the Skyrme model using the effective Lagrangian 
in flat space, R 3 , defined as 

C = f -f X 2 TrL,L» + ^-Tr\L,L v } 2 - ^d^X ~ jx\^ X 4 - 1), (3) 

which has the correct QCD scaling law. The glueball field x is normalized as 
x/xo — 1> C an d B are related to Xo an d the bag constant (B B = Bx 4 /4) re- 
spectively. Some years ago, Reinhardt and DangO] investigated the behav- 



ior of the skyrmion in the dense baryonic matter using the same Lagrangian 
on S 3 . They found that there is a phase transition into delocalized Skyrmion 
triggered by the vanishing gluon condensate. But, it is not easy to read out 
the order of the phase transition and the critical radius from their calcula- 
tions. In this work, we carry out the numerical analysis in more detail to 
investigate the characteristics of phase transition, particularly to see whether 
it is first or second order phase transition. 

We use the hedgehog ansatz for the skyrmion field configuration on S 3 , 

[/ = exp[fr-f(M)/G")], (4) 

where (/i, 8, <fi) are polar coordinates on S 3 (p) with < yU, 9 < it and < < 
27r and the metric is given by 

ds 2 = p 2 (d/j 2 + sin 2 /j,d8 2 + sin 2 /i sin 2 6d(j) 2 ). (5) 

The static energy of the system can be obtained from eq.(U) using hedgehog 
ansatz for chiral field and classical field for x(/i), 



E = ^J^dVC, 



2e, 

sin 2 / , 



^-{E 2 + E 4 + E x ), 



E 2 = 4vrL r d/isin 2 /iX 2 (7' 2 + 2 
Jo 
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• 2 /' 

sin fj, 
E x = Stt^LjT ^sinV{fx /2 H-L 2 [| X 4 (ln X 4 -l)H-|]}, (6) 

where (7 = C/i^esf-,,) 2 , B = I?/^/^) 4 and L is the dimensionless radius 
defined as L = e s fixP- 

Extremizing the static energy, we get the Euler equations for x(aO an d 
/(A*), 

/W + T^) + 2^/'(xx' + X 2 cot a*) 



L sin fx 



sin 2/ 2 1 sin / 



sin n L sin /x 

x " + 2cot/i X ' - ±{4L 2 B X 3 \nx+\x(f 2 + 2^)] = 0.. (8) 

c e; sin // 

where the chiral profile function and the glueball field satisfy the boundary 
conditions, for the B = 1 baryon, 

/(0) = 0, /(*)=*, x'(o) = o, xV) = o. (9) 

One can easily see that the solutions of the identity map /(//) = [J, with 
x(/i) = are the solutions of the above Euler equations for any radius, L. 
They are the delocalized solutions with no bag formation ||12| . 

We proceed to find numerical solutions with a typical parameter set [p~2f : 



e 4 = 4.5, U = 93MeV, C = 3.29 x !0 2 (MeV) 2 , B = (220MeV) 4 /e. (10) 

We can find two numerical solutions which are localized in addition to the 
identity map solution above mentioned. In Fig. 1, We plot the static energies 
of the solutions with respect to the radius of the hypersphere, L It is found 
that the two set of the localized solutions(curves (a) and (c) in Fig. 1) are 
possible for L > L t {= 1.9). One set of the solutions (curve (c) in Fig. 1) is 



not of our interest, Since we are interested in the lower energy configurations. 
At L = L c {= 1.9) the energy of the identity map solutions(curve(b)) is equal 
to the energy of the localized solutions (curve (a)). For a smaller size of the 
hypersphere, L < L c , the identity map solutions(delocalized phase) are the 
lower energy solutions while for the larger hypersphere, L > L c , the localized 
solutions are the minimum energy solutions. Hence there is a phase transi- 
tion from the localized to the delocalized phase at L c , which is of first order. 
The characteristic of the phase transition shown in Fig.l is quite different 
from the results without x field in the effective Lagrangian, where the phase 
transition is smooth and found to be of second order (L c = V2). Also the 
changes of the configurations of f(fi) and xif 1 ) with L are found out to be 
abrupt at L c as shown in Fig. 2(a) and 2(b) respectively. Since the delocal- 
ized solutions characterized by x(/ i ) = everywhere correspond to vanishing 
gluon condensate, the phase transition at L c from the localized phase with 
nonvanishing x{^) at /z = 7r to delocalized phase implies a discontinuous 
change of the gluon condensate at L c . These results can be considered as 



evidences of the first order phase transition WM triggered by the glueball field 



The order parameter, (a) , which is defined by 

(cr) = / dfi sin 2 fi cos f([i). (11) 

IX JO 

has been introduced as a useful quantity which effectively measures the chiral 
symmetry breaking. As shown in Fig. 3, one can see the discontinuous change 
of the order parameter (cr)o at L c . This suggests that the chiral symmetry is 
restored at the same radius L c and the phase transition is of first order. 

Some remarks are in order for the localized solutions. It is found that 
the localized solutions(the curves (a) and (c) in Fig.l) merge at L t and the 
localized solutions are no more possible for L < L t in numerical calculation. 
In order to see whether it is an numerical artifact, we look into the gradient 
of the profile function at /i = 0, f'{0) as a function of L. As shown in Fig. 
4, /'(0)'s of the two localized solutions are smoothly joined at L t and con- 
tinuously reverse their slopes respect to L, dL ■ This observation provides 
at least an partial explanation of the merging behavior of the two set of the 
localized solutions and the absence of localized solutions for L < L t other 
than the identity map solutions. 

It is also interesting to note that our analyses indicate that the metastable 
baryonic state, (a) in Fig. 1, may possibly exist between L t < L < L c , since 



the energy of this state is not much bigger than the identity map solution 
(b) in Fig. 1. 

In conclusion, in the frame work of the Skyrme model on S 3 with scale 
symmetry implemented, we find the evidence for the first order phase transi- 
tion at the critical radius of the hypersphere, L c ; from the localized phase to 
the delocalized phase, the chiral restoration and the gluon decondensation. 
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Figure Captions 

Fig. 1 The static energies with respect to the dimensionless radius, L for 
two set of the localized numerical solutions, (a) and (c), and for the identity 
map solutions, (b). 

Fig. 2 The numerical solutions of the chiral profile function /(/x), (a) and 
the glueball field x(/-0) (b)- 

Fig. 3 The order parameter (cr)o- The dahed line is for a model without 
glue ball field and the solid line for the extended model which is implemented 
by the glueball field. 

Fig. 4 The slops of the chiral profile function at \x — 0, /'(0) : The curve 
(a) corresponds to the lower energy localized solutions, the curve (a) in Fig. 
1, and the curve (c) corresponds to the higher energy localized solutions, the 
curve (c) in Fig. 1. 
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Fig. 1 
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Fig. 2 (a; 
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Fig. 2 (b) 
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Fig. 4 
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